Competition between Kondo and RKKY correlations in the presence of strong 
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We propose that competition between Kondo and magnetic correlations results in a novel univer- 
sality class for heavy fermion quantum criticality in the presence of strong randomness. Starting from 
an Anderson lattice model with disorder, we derive an effective local field theory in the dynamical 
mean-field theory (DMFT) approximation, where randomness is introduced into both hybridization 
and Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions. Performing the saddle-point analysis 
in the U(f ) slave-boson representation, we reveal its phase diagram which shows a quantum phase 
transition from a spin liquid state to a local Fermi liquid phase. In contrast with the clean limit 
of the Anderson lattice model, the effective hybridization given by holon condensation turns out 
to vanish, resulting from the zero mean value of the hybridization coupling constant. However, we 
show that the holon density becomes finite when variance of hybridization is sufficiently larger than 
that of the RKKY coupling, giving rise to the Kondo effect. On the other hand, when the variance 
of hybridization becomes smaller than that of the RKKY coupling, the Kondo effect disappears, re- 
sulting in a fully symmetric paramagnetic state, adiabatically connected with the spin liquid state of 
the disordered Ifeisenberg model. We investigate the quantum critical point beyond the mean-field 
approximation. Introducing quantum corrections fully self-consistently in the non-crossing approx- 
imation, we prove that the local charge susceptibility has exactly the same critical exponent as the 
local spin susceptibility, suggesting an enhanced symmetry at the local quantum critical point. This 
leads us to propose novel duality between the Kondo singlet phase and the critical local moment 
state beyond the Landau-Ginzburg- Wilson paradigm. The Landau-Ginzburg- Wilson forbidden du- 
ality serves the mechanism of electron fractionalization in critical impurity dynamics, where such 
fractionalized excitations are identified with topological excitations. 

PACS numbers: 



I. INTRODUCTION 

Interplay between interactions and disorders has been 
one of the central issues in modern condensed matter 
physics [H, [3 ■ In the weakly disordered metal the lowest- 
order interaction-correction was shown to modify the 
density of states at the Fermi energy in the diffusive 
regime giving rise to non-Fermi liquid physics par- 
ticulary in low dimensions less than d — Z while fur- 
ther enhancement of electron correlations was predicted 
to cause ferromagnetism In an insulating phase spin 
glass appears ubiquitously, where the average of the spin 
moment vanishes in the long time scale, but local spin 
correlations become finite, making the system away from 
equilibrium [5]. 

An outstanding question is the role of disorder in the 
vicinity of quantum phase transitions [1, 0] , where effec- 
tive long-range interactions associated with critical fluc- 
tuations appear to cause no n- Fermi liquid physics 0, Q . 
Unfortunately, complexity of this problem did not allow 
comprehensive understanding until now. In the vicinity 
of the weakly disordered ferromagnetic quantum critical 
point, an electrical transport-coefficient has been stud- 
ied, where the crossover temperature from the ballistic 
to diffusive regimes is much lowered due to critical fluc- 
tuations, compared with the disordered Fermi liquid 9]. 
Generally speaking, the stability of the quantum critical 
point should be addressed, given by the Harris criterion 



[lO|. When the Harris criterion is not satisfied, three 
possibilities are expected to arise [7]. The first two pos- 
sibilities are emergence of new fixed points, associated 
with either a finite-randomness fixed point satisfying the 
Harris criterion at this new fixed point or an infinite ran- 
domness fixed point exhibiting activated scaling behav- 
iors. The last possibility is that quantum criticality can 
be destroyed, replaced with a smooth crossover. In addi- 
tion, even away from the quantum critical point the disor- 
dered system may show non-universal power law physics, 
called the Griffiths phase [ll| . Effects of rare regions are 
expected to be strong near the inflnite randomness fixed 
point and the disorder-driven crossover region Q . 

This study focuses on the role of strong randomness 
in the heavy fermion quantum transition. Heavy fermion 
quantum criticality is believed to result from competition 
between Kondo and RKKY (Ruderman-Kittel-Kasuya- 
Yosida) interactions, where larger Kondo couplings give 
rise to a heavy fermion Fermi liquid while larger RKKY 
interactions cause an antiferromagnetic metal 0, [1, [13 ■ 
Generally speaking, there are two competing view points 
for this problem. The first direction is to regard the 
heavy fermion transition as an antiferromagnetic transi- 
tion, where critical spin fluctuations appear from heavy 
fermions. The second view point is that the transi- 
tion is identified with breakdown of the Kondo effect, 
where Fermi surface fluctuations are critical excitations. 
The flrst scenario is described by the Hertz-Moriya-Millis 
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(HMM) theory in terms of heavy electrons coupled with 
antiferromagnetic spin fluctuations, the standard model 
for quantum criticality [l^ . There are two ways to realize 
the second scenario depending on how to describe Fermi 
surface fluctuations. The first way is to express Fermi 
surface fluctuations in terms of a hybridization order pa- 
rameter called holon in the slave-boson context [13, [3] ■ 
This is usually referred as the Kondo breakdown scenario. 
The second one is to map the lattice problem into the 
single site one resorting to the dynamical mean-field the- 
ory (DMFT) approximation (l6j , where order parameter 
fluctuations are critical only in the time direction. This 
description is called the locally critical scenario [l7j . 

Each scenario predicts its own critical physics. Both 
the HMM theory and the Kondo breakdown model are 
based on the standard picture that quantum criticality 
arises from long-wave-length critical fluctuations while 
the locally quantum critical scenario has its special struc- 
ture, that is, locally (space) critical (time). Critical fluc- 
tuations are described by z = 2 in the HMM theory due 
to finite-wave vector ordering while by z = 3 in the 
Kondo breakdown scenario associated with uniform "or- 
dering" [l^ , where z is the dynamical exponent express- 
ing the dispersion relation for critical excitations. Thus, 
quantum critical physics characterized by scaling expo- 
nents is completely different between these two models. 
In addition to qualitative agreements with experiments 
depending on compounds , these two theories do not al- 
low the ijj /T scaling in the dynamic susceptibility of their 
critical modes because both theories live above their up- 
per critical dimensions. On the other hand, the locally 
critical scenario gives rise to the a;/T scaling behavior 
for the dynamic spin susceptibility |I7l | while it seems to 
have some difficulties associated with some predictions 
for transport coefRcients. 

We start to discuss an Ising model with Gaussian ran- 
domness for its exchange coupling, called the Edwards- 
Anderson model • Using the replica trick and perform- 
ing the saddle-point analysis, one can find a spin glass 
phase when the average value of the exchange interaction 
vanishes, characterized by the Edwards- Anderson order 
parameter without magnetization. Applying this con- 
cept to the Heisenberg model with Gaussian randomness, 
quantum fluctuations should be incorporated to take into 
account the Berry phase contribution carefully. It was 
demonstrated that quantum corrections in the DMFT 
approximation lead the spin glass phase unstable at fi- 
nite temperatures, resulting in a spin liquid state when 
the average value of the exchange coupling vanishes [l^ ■ 
It should be noted that this spin liquid state differs from 
the spin liquid phase in frustrated spin systems in the re- 
spect that the former state originates from critical single- 
impurity dynamics while the latter phase results from 
non-trivial s pat ial spin correlations described by gauge 
fluctuations [19| . The spin liquid phase driven by strong 
randomness is characterized by its critical spin spectrum, 
given by the uj/T scaling local spin susceptibility [T8| . 

Introducing hole doping into the spin liquid state. Par- 



collet and Georges examined the disordered t-J model 
within the DMFT approximation Using the U(f) 

slave-boson representation, they found marginal Fermi- 
liquid phenomenology, where the electrical transport is 
described by the T-linear resistivity, resulting from the 
marginal Fermi-liquid spectrum for collective modes, 
here the uj/T scaling in the local spin susceptibility. 
They tried to connect this result with physics of high 
Tc cuprates. 

In this study we introduce random hybridization with 
conduction electrons into the spin liquid state. Our orig- 
inal motivation was to explain both the oj/T scaling in 
the smn spectrum [ST] and the typical T-linear resistiv- 
ity |22| near the heavy fcrmion quantum critical point. 
In particular, the presence of disorder leads us to the 
DMFT approximation naturally [23|], expected to result 
in the u/T scaling for the spin spectrum [Tsj . 

Starting from an Anderson lattice model with disorder, 
we derive an effective local field theory in the DMFT ap- 
proximation, where randomness is introduced into both 
hybridization and RKKY interactions. Performing the 
saddle-point analysis in the U(I) slave-boson representa- 
tion, we reveal its phase diagram which shows a quantum 
phase transition from a spin liquid state to a local Fermi 
liquid phase. In contrast with the clean limit of the An- 
derson lattice model [f3, HI] , the effective hybridization 
given by holon condensation turns out to vanish, result- 
ing from the zero mean value of the hybridization cou- 
pling constant. However, we show that the holon density 
becomes finite when variance of hybridization is suffi- 
ciently larger than that of the RKKY coupling, giving 
rise to the Kondo effect. On the other hand, when the 
variance of hybridization becomes smaller than that of 
the RKKY coupling, the Kondo effect disappears, result- 
ing in a fully symmetric paramagnetic state, adiabatically 
connected with the spin liquid state of the disordered 
Heisenberg model [l8| . 

Our contribution compared with the previous works 
[2^ 1 is to introduce RKKY interactions between localized 
spins and to observe the quantum phase transition in the 
heavy fermion system with strong randomness. The pre- 
vious works focused on how the non-Fermi liquid physics 
can appear in the Kondo singlet phase away from quan- 
tum criticality [IJ]. A huge distribution of the Kondo 
temperature Tk turns out to cause such non- Fermi liquid 
physics, originating from the finite density of unscreened 
local moments with almost vanishing Tk, where the Tk 
distribution may result from either the Kondo disorder 
for localized electrons or the proximity of the Anderson 
localization for conduction electrons. Because RKKY in- 
teractions are not introduced in these studies, there al- 
ways exist finite Tk contributions. On the other hand, 
the presence of RKKY interactions gives rise to break- 
down of the Kondo effect, making Tk — identically in 
the strong RKKY coupling phase. 

In Ref. [2^ the role of random RKKY interactions was 
examined, where the Kondo coupling is fixed while the 
chemical potential for conduction electrons is introduced 
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as a random variable with its variance W. Increasing the 
randomness of the electron chemical potential, the Fermi 
liquid state in W < Wc turns into the spin liquid phase 
in > Wc, which displays the marginal Fermi- liquid 
phenomenology due to random RKKY interactions 25], 
where the Kondo effect is suppressed due to the proxim- 
ity of the Anderson localization for conduction electrons 
[24| . However, the presence of finite Kondo couplings still 
gives rise to Kondo screening although the Tk distribu- 
tion differs from that in the Fermi liquid state, associated 
with the presence of random RKKY interactions. In ad- 
dition, the spin liquid state was argued to be unstable 
against the spin glass phase at low temperatures, maybe 
resulting from the fixed Kondo interaction. On the other 
hand, we do not take into account the Anderson local- 
ization for conduction electrons, and introduce random 
hybridization couplings. As a result, the Kondo effect 
is completely destroyed in the spin liquid phase, thus 
quantum critical physics differs from the previous study 
of Ref. . In addition, the spin liquid phase is stable 
at finite temperatures in the present study [T8l|. 

We investigate the quantum critical point beyond the 
mean-field approximation. Introducing quantum correc- 
tions fully self-consistently in the non-crossing approxi- 
mation [26| , we prove that the local charge susceptibility 
has exactly the same critical exponent as the local spin 
susceptibility. This is quite unusual because these cor- 
relation functions are symmetry-unrelated in the lattice 
scale. This reminds us of deconfined quantum criticality 
[27| . where the Landau-Ginzburg- Wilson forbidden con- 
tinuous transition may appear with an enhanced emer- 
gent symmetry. Actually, the continuous quantum tran- 
sition was proposed between the antiferromagnetic phase 
and the valence bond solid state [27[. In the vicinity of 



the quantum critical point the spin-spin correlation func- 
tion of the antiferromagnetic channel has the same scal- 
ing exponent as the valence-bond correlation function, 
suggesting an emergent 0(5) symmetry beyond the sym- 
metry 0(3) XZ4 of the lattice model [28] and confirmed 
by the Monte Carlo simulation of the extended Heisen- 
berg model [2^. Tanaka and Hu proposed an effective 
0(5) nonlinear a model with the Wess-Zumino-Witten 
term as an effective field theory for the Landau- Ginzburg- 
Wilson forbidden quantum critical point [l^ , expected to 
allow fractionalized spin excitations due to the topologi- 
cal term. This proposal can be considered as generaliza- 
tion of an antiferromagnetic spin chain, where an effective 
field theory is given by an 0(4) nonlinear a model with 
the Wess-Zumino-Witten term, which gives rise to frac- 
tionalized spin excitations called spinous, identified with 
topological solitons [s^]. Applying this concept to the 
present quantum critical point, the enhanced emergent 
symmetry between charge (holon) and spin (spinous) lo- 
cal modes leads us to propose novel duality between the 
Kondo singlet phase and the critical local moment state 
beyond the Landau-Ginzburg- Wilson paradigm. We sug- 
gest an 0(4) nonlinear a model in a nontrivial manifold 
as an effective field theory for this local quantum critical 



point, where the local spin and charge densities form an 
0(4) vector with a constraint. The symmetry enhance- 
ment serves the mechanism of electron fractionalization 
in critical impurity dynamics, where such fractionalized 
excitations are identified with topological excitations. 

This paper is organized as follows. In section II we 
introduce an effective disordered Anderson lattice model 
and perform the DMFT approximation with the replica 
trick. Equation (j4]) is the main result in this section. In 
section III we perform the saddle-point analysis based on 
the slave-boson representation and obtain the phase dia- 
gram showing breakdown of the Kondo effect driven by 
the RKKY interaction. We show spectral functions, self- 
energies, and local spin susceptibility in the Kondo phase. 
Figures (l)-(3) with Eqs. ([IH1)-(I1I]) and (I121)-(IM1) are 
main results in this section. In section IV we investigate 
the nature of the impurity quantum critical point based 
on the non-crossing approximation beyond the previous 
mean-field analysis. We solve self-consistent equations 
analytically and find power-law scaling solutions. As a 
result, we uncover the marginal Fermi-liquid spectrum 
for the local spin susceptibility. We propose an effective 
field theory for the quantum critical point and discuss 
the possible relationship with the deconfined quantum 
critical point. In section V we summarize our results. 

The present study extends our recent publication [sj, 
showing both physical and mathematical details. 



II. AN EFFECTIVE DMFT ACTION FROM AN 
ANDERSON LATTICE MODEL WITH STRONG 
RANDOMNESS 

We start from an effective Anderson lattice model 

ij,a i<7 

+ J.,S, • S, + Y.^v4j,, + H.c), (1) 



electrons and 



is a hopping integral for conduction 
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are random RKKY and hybridization coupling constants, 
respectively. Here, M is the spin degeneracy and z is 
the coordination number. Randomness is given by the 
Gaussian distribution 



e7=0. 
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The disorder average can be performed in the replica 
trick . Performing the disorder average in the Gaussian 
distribution function, we reach the following expression 
for the replicated effective action 
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where tr, ct' = 1, M is the spm index and a,b — 1, n 
is the rephca index. In appendix A we derive this reph- 
cated action from Eq. ([1]). 

One may ask the role of randomness of Ed, generating 

- f^dr f dr' J2 4"(r)<.(r)4^,(r')dt(r'), 

■JO JO /J, 
iaa'ab 

where density fluctuations are involved. This contribu- 
tion is expected to support the Kondo effect because such 
local density fluctuations help hybridization with conduc- 
tion electrons. In this paper we fix Ed as a constant value 
in the Kondo limit, allowed as long as its variance is not 
too large to overcome the Kondo limit. 

One can introduce randomness in the hopping integral 
of conduction electrons. But, this contribution gives rise 
to the same effect as the DMFT approximation in the 
z — 7> CXI Bethe lattice [l^]. In this respect randomness 
in the hopping integral is naturally introduced into the 
present DMFT study. 

The last disorder contribution can arise from random- 
ness in the electron chemical potential, expected to cause 
the Anderson localization for conduction electrons. Actu- 
ally, this results in the metal-insulator transition at the 
critical disorder strength, suppressing the Kondo effect 
in the insulating phase. Previously, the Griffiths phase 
for non-Fermi liquid physics has been attributed to the 
proximity effect of the Anderson localization [2J| . In this 
work we do not consider the Anderson localization for 
conduction electrons. 

We observe that the disorder average neutralizes spa- 
tial correlations except for the hopping term of conduc- 
tion electrons. This leads us to the DMFT formulation, 
resulting in an effective local action for the strong ran- 
dom Anderson lattice model 
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where G'^\j^^{r — r') is the local Green's function for 

conduction electrons and Rp^'^jsi''' "''"') the local spin 
susceptibility for localized spins, given by 



fiab I _ , 
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= -(r.[C(r)4^(r')]), 



(5) 



respectively. Eq. Q with Eq. ([5]) serves a completely 
self-consistent framework for this problem. Derivation of 
Eq. ^ from Eq. ^ is shown in appendix B. 

This effective model has two well known lim its, cor- 
responding to the disordered Heisenberg model and 
the disordered Anderson lattice model without RKKY 
interactions [l^], respectively. In the former case a spin 
liquid state emerges due to strong quantum fluctuations 
while a local Fermi liquid phase appears at low temper- 
atures in the latter case as long as the Tk distribution 
is not so broadened enough. In this respect it is natural 
to consider a quantum phase transition driven by the ra- 
tio between variances for the RKKY and hybridization 
couplings. 



III. PHASE DIAGRAM 

A. Slave boson representation and mean field 
approximation 

We solve the effective DMFT action based on the U(I) 
slave boson representation 



dZ 



Sla' = fa^fa' ~ lo^a 



(6) 
(7) 



with the single occupancy constraint -I- 
fa{r)f^{T) = 1, where = Y.a fa^falM. 
In the mean field approximation we replace the holon 
operator b°' with its expectation value (6°) = 6". Then, 
the effective action Eq. ^ becomes 
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where A" is a lagrange multiplier field to impose the con- ZZr* '^^ ( \h'^h'' (^')\ 

straint and = (/IVa). ^ ^c/ <r<r' l^jo o , U^J 

Taking the Af — ?> cxi limit, we obtain self-consistent 
equations for self-energy corrections, 

respectively, where local Green's functions are given by 
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1/2 In the paramagnetic and symmetric replica phase these 

~j^'^fc(7a'\~)\" "7 ' ^-"^7 Green's functions are diagonal in the spin and replica in- 

2 dices, i.e., G°^„,(t) = (5o&^<T<T'Ga;(T) with x = c, f,cf,fc. 
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in the frequency space. Note that n is the replica index 



V'^ 2 ^'^'^ 1^^^ terms in Eqs. (f20| -(|2T |) vanish in the limit of 
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The constraint equation is given by 



\b\'+j2{nu)^i. 



(24) 



The main difference between the clean and disordered 
cases is that the off diagonal Green's function GfdioJi) 
should vanish in the presence of randomness in V with 
its zero mean value while it is proportional to the con- 
densation b when the average value of V is finite. In 
the present situation we find &° — (f^'^c^) — while 
ib'^yb" = {f^^clc^^Jl.) = \b\Hai, ^0. As a resuh, Eqs. 
([20|l and ([21]) are identically vanishing in both left and 
right hand sides. This implies that the Kondo phase 
is not characterized by the holon condensation but de- 
scribed by finite density of holons. It is important to no- 
tice that this gauge invariant order parameter does not 
cause any kinds of symmetry breaking for the Kondo ef- 
fect as it should be. 



B. Numerical analysis 

We use an iteration method in order to solve the mean 
field equations ([111), ([111), ([201), ([HI), (113, and ([Ml For 
a given E^i ~\- A, we use iterations to find all Green's func- 
tions from Eqs. (ITH1)-(11T|) with Eq. and 6^ from 
Eq. Then, we use Eq. ^ to calculate A and Ed- 
We adjust the value of i^jj -I- A in order to obtain the de- 
sirable value for Ed- Using the obtained A and b^, we 
calculate the Green's functions in the real frequency by 
iterations. In the real frequency calculation we introduce 
the following functions [32] 



a±{t) = / duje-'^'pf{Lu)f{±LolT), 



(25) 



where Pfito) — —lniGf{uj + iO^)/7T is the density of states 
for f-electrons, and f{x) — l/(exp(x) -I- 1) is the Fermi- 
Dirac distribution function. Then, the self-energy cor- 
rection from spin correlations is expressed as follows 



Sj(za;;) 



2M 



Jo 
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Performing the Fourier transformation, we calculate 
a±{t) and obtain Ej(a;). 

Figure [T] shows the phase diagram of the strongly dis- 
ordered Anderson lattice model in the plane of {V,J), 
where V and J are variances for the Kondo and RKKY 




FIG. 1: The phase diagram of the strongly disordered An- 
derson lattice model in the DMFT approximation {Ed = — 1, 
11 = 0, r = 0.01, t= 1, M = 2). 



interactions, respectively. The phase boundary is charac- 
terized by l&p =0, below which |6p ^0 appears to cause 
effective hybridization between conduction electrons and 
localized fermions although our numerical analysis shows 
if^Ccr) = 0, meaning i;c/(/c)(«w) = and Gcf{fc){i^^) = 



in Eqs. ^ and ([21]) 

In Fig. [H one finds that the effective hybridization 
enhances the scattering rate of conduction electrons dra- 
matically around the Fermi energy while the scattering 
rate for localized electrons becomes reduced at the reso- 
nance energy. Enhancement of the imaginary part of the 
conduction-electron self-energy results from the Kondo 
effect. In the clean situation it is given by the delta 
function associated with the Kondo effect [2g . This self- 
energy effect refiects the spectral function, shown in Fig. 
|3l where the pseudogap feature arises in conduction elec- 
trons while the sharply defined peak appears in localized 
electrons, identified with the Kondo resonance although 
the description of the Kondo effect differs from the clean 
case. Increasing the RKKY coupling, the Kondo effect 
is suppressed as expected. In this Kondo phase the local 
spin susceptibility is given by Fig. 21 displaying the typ- 
ical oj-linear behavior in the low frequency limit, nothing 
but the Fermi liquid physics for spin correlations [20| . 
Increasing J, incoherent spin correlations are enhanced, 
consistent with spin liquid physics [20| . 

One can check our calculation, considering the J = 
limit to recover the known result. In this limit we obtain 
an analytic expression for Vc at half filling (/i = 0) 



Ed 



(27) 



Pr = 



dujpo(u}) , (28) 



where po{i^) — -^Vl—oj^ is the bare density of states of 
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FIG. 2; The imaginary part of the self-energy of conduction FIG. 3: Density of states of conduction (pc(w)) and localized 
electrons and that of localized electrons for various values of (P/(^)) electrons for various values oi J {V = 0.5, Ed = -0.7, 
J {V = 0.5, Ea = -0.7, = 0, T = 0.01, t = 1, M=2). ^ = 0, T = 0.01, t = l, M = 2). 



conduction electrons. One can check Vc{J = 0) — >■ in 
the zero temperature limit because Pc — ^ oo. 

IV. NATURE OF QUANTUM CRITICALITY 

A. Beyond the saddle-point analysis : Non-crossing 
approximation 

Resorting to the slave-boson mean-field approxima- 
tion, wc discussed the phase diagram of the strongly dis- 
ordered Anderson lattice model, where a quantum phase 
transition appears from a spin liquid state to a dirty 
" heavy- fcrmion" Fermi liquid phase, increasing V/J, the 
ratio of variances of the hybridization and RKKY in- 
teractions. Differentiated from the heavy-fermion quan- 
tum transition in the clean situation, the order parameter 
turns out to be the density of holons instead of the holon 
condensation. Evaluating self-energies for both conduc- 
tion electrons and localized electrons, we could identify 
the Kondo effect from each spectral function. In addi- 
tion, wc obtained the local spin susceptibility consistent 
with the Fermi liquid physics. 

The next task will be on the nature of quantum crit- 
icality between the Kondo and spin liquid phases. This 
question should be addressed beyond the saddle-point 



analysis. Introducing quantum corrections in the non- 
crossing approximation, justified in the M — > oo limit, 
we investigate the quantum critical point, where density 
fluctuations of holons are critical. 

Releasing the slave-boson mean-field approximation to 
take into account holon excitations, we reach the follow- 
ing self-consistent equations for self-energy corrections, 

+ GrUrmta'Ar) + K's'sA-r)], 

ss' 

(30) 

sc '' 

(31) 

sc 

(32) 
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FIG. 4: Local spin susceptibility for various values of J (V = 
0.5, Ed = -0.7, fi^O,T = 0.01, t = 1, M = 2). 



^b\r) - ^ E ir)G,'^.'A-r). (33) 



Since we considered the paramagnetic and replica sym- 
metric phase, it is natural to assume such symmetries at 
the quantum critical point. Note that the off diagonal 
self-energies, I]c/(«w;) and I]/c(ia;;), are just constants 
and proportional to {.f^Ca) and {clfa), respectively. As 
a result, Y,cf{iwi) = Sjc(iw;) = should be satisfied at 
the quantum critical point as the Kondo phase because 
of {flca) = (cJ./(t) = 0. Then, we reach the following 
self-consistent equations called the non-crossing approx- 
imation 



M 



(34) 



S/(r) = —G,{T)G,{T)~J'[Gf{T)YGf{-T),{i^) 



Sb(r) = V^GrX-T)Gf{T). 
Local Green's functions are given by 



Gc{iuJi) 
Gfiiuji) 
Gb{m) 



iuji+ fi~ T,ciiuji) 



-1 



lUJl - 



(36) 

(37) 
(38) 
(39) 



where uji — {21 + 1)t:T is for fermions and vi = 21ttT is 
for bosons. 



B. Asymptotic behavior at zero temperature 

Calling quantum criticality, power-law scaling solu- 
tions are expected. Actually, if the second term is ne- 
glected in Eq. Eqs. (1551) and (1551) are reduced 



to those of the multi-channel Kondo effect in the non- 
crossing approximation ^26i] . Power-law solutions are well 
known in the regime of 1/Tk t /3 — 1/T ^ oo, 
where Tk = D[Tc/nD]^/'^'^ explnEd/MTc] is an effective 
Kondo temperature [33] with the conduction bandwidth 

D and effective hybridization Fc = npc-jj- In the pres- 
ence of the RKKY interaction [the second term in Eq. 
([55]) ]. the effective hybridization will be reduced, where 

Tc is replaced with F;,' irpdj^ — J^)- 
Our power-law ansatz is as follows 



Gc = 
Gf = 
Gb — 



(40) 
(41) 
(42) 



where Ac, Af, and Ab are positive numerical constants. 
In the frequency space these are 



Gfiio) = AfCAf-iUJ^'-\ 
Gb{uj) - AbGA,-iUj'^'-\ 

dx- 



(43) 
(44) 
(45) 



where Ca,,,,, - u^^^^^^. 

Inserting Eqs. (|131)-(|lSl) into Eqs. (gUl-dMl), we obtain 
scaling exponents of Ac, A/, and A^. In appendix C-1 
we show how to find such critical exponents in a detail. 
Two fixed points are allowed. One coincides with the 
multi-channel Kondo effect, given by Ac — 1, and — 



M 



Ab = 



with M = 2, where contributions from 



A/+1' "0 M+1 

spin fluctuations to self-energy corrections are irrelevant, 
compared with holon fluctuations. The other is Ac = 1 
and A/ = Ah ~ ^ , where spin correlations are critical as 
much as holon fluctuations. 

One can understand the critical exponent A/ = 1/2 as 
the proximity of the spin liquid physics [l8l | . Considering 
the V limit, we obtain the scaling exponents of 
Ac = 1 and Af = 1/2 from the scaling equations (jC22p 
and ((C23)) . Thus, Gc(a;) - sgn(cj) and Gf{uj) ^ 1/^ 
result for w — >■ 0. In this respect both spin fluctuations 
and holon excitations are critical as equal strength at this 
quantum critical point. 



C. Finite temperature scaling behavior 

We solve Eqs. dH!)-® in the regime r, /3 > 1/Tk 
with arbitrary t//3, where the scaling ansatz at zero tem- 
perature is generalized as follows 



Gc(r) = ^c/3-^^9c(^), 
Gfir) = Afr^fgf[^), 
Gbir) = Abr^^gb[j). 



(46) 
(47) 
(48) 
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sin(7ra::) 



(49) 



with a = c, /, 6 is the scaling function at finite tempera- 
tures. In the frequency space we obtain 





= A,p^- 




(50) 


Gf{iuji) 






(51) 


Gb{m) 


= A,p^- 




(52) 



where uji = {21 + l)7r, vi = 2171, and 
Jo 



(53) 



Inserting Eqs. dSHll-dSl]) into Eqs. dSU-dSS]), we find 
two fixed points, essentially the same as the case of T = 0. 
But, scaling functions of $c(*'^;), ^/(iwi), and ^^(jw;) 
are somewhat complicated. All scaling functions are de- 
rived in appendix C-2. 



D. Spin susceptibility 

We evaluate the local spin susceptibility, given by 



sin(7rr//3) 



2Af 



(54) 



The imaginary part of the spin susceptibility x (cj) = 
Im x(a; -h zO+) can be found from 



X(r) 



dw 



TT 1 - e-'3" 
Inserting the scaling ansatz 

into Eq. ([55]) with Eq. ([5l)) . we obtain 



(55) 



(56) 



This coincides with the spin spectrum of the spin liquid 
state when F = [13]. 

E. Discussion : Deconfined local quantum 
criticality 

The local quantum critical point characterized by Ac = 
1 and = Aj, = 1/2 is the genuine critical point in the 
spin-liquid to local Fermi-liquid transition because such 
a fixed point can be connected to the spin liquid state 
(Ac = 1 and A/ = 1/2) naturally. This fixed point re- 
sults from the fact that the spinon self-energy correction 
from RKKY spin fluctuations is exactly the same order as 
that from critical holon excitations. It is straightforward 
to see that the critical exponent of the local spin sus- 
ceptibility is exactly the same as that of the local charge 
susceptibility (2 Ay = 2 A;, = 1), proportional to l/r. 
Since the spinon spin-density operator differs from the 
holon charge-density operator in the respect of symme- 
try at the lattice scale, the same critical exponent implies 
enhancement of the original symmetry at low energies. 
The symmetry enhancement sometimes allows a topo- 
logical term, which assigns a nontrivial quantum num- 
ber to a topological soliton, identifled with an excitation 
of quantum number fractionalization. This mathemati- 
cal structure is actually realized in an antiferromagnetic 
spin chain [30| , generalized into the two dimensional case 

We propose the following local field theory in terms of 
physically observable fields 



where 



i:i*"(r)5 |*°(t)| 



9 
'2M 

top 7 



1 e- 



iff 



dT dT'*''^(T)Y''''(r-r')*''(T') 
Jo Jo 

(60) 



(61) 



2Af 



TT 1 — e ^ ^ ^ ' \sa\{'KT / 13) 
Changing the variable t = i(T//3 — 1/2), we obtain 



dx 



e^-^ \cosh(7rt) 
As a result, we find the scaling function 



2Af 



(f>{x) = 2(27r) 



represents an 0(4) vector, satisfying the constraint of the 
delta function. T''''(t — t') determines dynamics of the 
0(4) vector, resulting from spin and holon dynamics in 
principle. However, it is extremely difficult to derive Eq. 
([SO)) from Eq. @ because the density part for the holon 
field in Eq. ([60)) cannot result from Eq. ([4]) in a stan- 
dard way. What we have shown is that the renormalized 
dynamics for the 0(4) vector field follows l/r asymptot- 
ically, where r is the imaginary time. This information 
should be introduced in T''''(r — r'). g (x V/ J is an effec- 
2Af-i ^^^Yi (—\ -*^(^/ ix/2'K)T{Af — ix/2Ti) ^jyg coupling constant, and Stop is a possible topological 
.2/ r(2A/) ■ term. 

(59) One can represent the 0(4) vector generally as follows 



(57) 



(58) 
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*° : T — > (sin 9" (t) sin (P^ir) cos ip^ir), sin 9" (t) 



where 9'^{t), 0°(t), (P°'{t) are three angle coordinates for 
the 0(4) vector. It is essential to observe that the target 
manifold for the 0(4) vector is not a simple sphere type, 
but more complicated because the last component of the 
0(4) vector is the charge density field, where three spin 
components lie in -1 < S^{t), S^{t), S^{t) < 1 while 
the charge density should be positive, < p°'{t) < 1. 
This leads us to identify the lower half sphere with the 
upper half sphere. Considering that sin0"(T) can be 
folded on 7r/2, we are allowed to construct our target 
manifold to have a periodicity, given by ^"^{9"', ip°',Lp°-) — 
^"(tt — 0", (j)"', If"-). This folded space allows a nontrivial 
topological excitation. 

Suppose the boundary configuration of 
*"(O,0°,V3'';t = 0) and *'^(7r, <?!)'^, r = ^), con- 
nected by *"(7r/2,(^",93'';0 < r < /3). Interestingly, 
this configuration is topologically distinguishable 
from the configuration of ^°(0, 0", 1^9"; r = 0) and 
*''(O,0'', r = P) with *'^(7r/2, 0"^, (/j-^; < r < /3) 
because of the folded structure. The second configura- 
tion shrinks to a point while the first excitation cannot, 
identified with a topologically nontrivial excitation. This 
topological excitation carries a spin quantum number 
1/2 in its core, given by ^°-{n/2, , ip°-]0 < r < /3) = 

^sin 0" (r ) cos v?" (t) , sin 0" (r) sin </'"(''"), cos i/)'^ (r) , 0^ 

This is the spinon excitation, described by an 0(3) 
nonlinear a model with the nontrivial spin correlation 
function Y°-^{t — r'), where the topological term is 
reduced to the single spin Berry phase term in the 
instanton core. 

In this local impurity picture the local Fermi liquid 
phase is described by gapping of instantons while the spin 
liquid state is characterized by condensation of instan- 
tons. Of course, the low dimensionality does not allow 
condensation, resulting in critical dynamics for spinous. 
This scenario clarifies the Landau-Ginzburg- Wilson for- 
bidden duality between the Kondo singlet and the critical 
local moment for the impurity state, allowed by the pres- 
ence of the topological term. 

If the symmetry enhancement does not occur, the ef- 
fective local field theory will be given by 

+ ^S^ir)R^\r - T')S'{r')] + Sb (63) 



_ (f)" (r) sin ifi" (r) , sin 9" (r) cos (j)" (r) , cos 9" (r)) , (62) 
I 

with the single-spin Berry phase term 

SB^-2mS du dT—S''{u,T)duS"{u,T)xdrS"{u,T), 
Jo Jo 47r 

where charge dynamics x''''(t — t') will be different from 
spin dynamics R°'^{t — t'). This will not allow the 
spin fractionalization for the critical impurity dynamics, 
where the instanton construction is not realized due to 
the absence of the symmetry enhancement. 



V. SUMMARY 

In this paper we have studied the Anderson lattice 
model with strong randomness in both hybridization and 
RKKY interactions, where their average values are zero. 
In the absence of random hybridization quantum fluctua- 
tions in spin dynamics cause the spin glass phase unstable 
at finite temperatures, giving rise to the spin liquid state, 
characterized by the lu/T scaling spin spectrum consis- 
tent with the marginal Fermi- liquid phenomenology pLSj . 
In the absence of random RKKY interactions the Kondo 
effect arises , but differentiated from that in the clean 
case. The dirty "heavy fermion" phase in the strongly 
disordered Kondo coupling is characterized by a finite 
density of holons instead of the holon condensation. But, 
effective hybridization exists indeed, causing the Kondo 
resonance peak in the spectral function. As long as vari- 
ation of the effective Kondo temperature is not so large, 
this disordered Kondo phase is identified with the local 
Fermi liquid state because essential physics results from 
single impurity dynamics, differentiated from the clean 
lattice model. 

Taking into account both random hybridization and 
RKKY interactions, we find the quantum phase transi- 
tion from the spin liquid state to the local Fermi liquid 
phase at the critical {Vc, Jc)- Each phase turns out to be 
adiabatically connected with each limit, i.e., the spin liq- 
uid phase when V = and the local Fermi liquid phase 
when J — 0, respectively. Actually, we have checked this 
physics, considering the local spin susceptibility and the 
spectral function for localized electrons. 

In order to investigate quantum critical physics, we 
introduce quantum corrections from critical holon fluc- 
tuations in the non-crossing approximation beyond the 
slave-boson mean-fleld analysis. We find two kinds of 
power-law scaling solutions for self-energy corrections of 
conduction electrons, spinous, and holons. The first solu- 
tion turns out to coincide with that of the multi-channel 
Kondo effect, where effects of spin fiuctuations are sub- 
leading, compared with critical holon fluctuations. In 
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this respect this quantum critical point is characterized 
by breakdown of the Kondo effect while spin fluctuations 
can be neglected. On the other hand, the second scaling 
solution shows that both holon excitations and spinon 
fluctuations are critical as the same strength, reflected in 
the fact that the density-density correlation function of 
holons has the exactly the same critical exponent as the 
local spin-spin correlation function of spinous. 

We argued that the second quantum critical point im- 
plies an enhanced emergent symmetry from 0(3) x 0(2) 
(spin(X)charge) to 0(4) at low energies, forcing us to con- 
struct an 0(4) nonlinear a model on the folded target 
manifold as an effective field theory for this disorder- 
driven local quantum critical point. Our effective local 
field theory identifies spinous with instantons, describ- 
ing the local Fermi-liquid to spin-liquid transition as the 
condensation transition of instantons although dynamics 
of instantons remains critical in the spin liquid state in- 
stead of condensation due to low dimensionality. This 
construction completes novel duality between the Kondo 
and critical local moment phases in the strongly disor- 
dered Anderson lattice model. 

We explicitly checked that the similar result can be 
found in the extended DMFT for the clean Kondo lat- 
tice model, where two fixed point solutions are allowed 
[3^ Issj . One is the same as the multi-channel Kondo 
effect and the other is essentially the same as the second 
solution in this paper. In this respect we believe that the 
present scenario works in the extended DMFT framework 
although applicable to only two spatial dimensions [T7| . 

One may suspect the applicability of the DMFT frame- 
work for this disorder problem. However, the hybridiza- 
tion term turns out to be exactly local in the case of 
strong randomness while the RKKY term is safely ap- 
proximated to be local for the spin liquid state, expected 
to be stable against the spin glass phase in the case of 
quantum spins. This situation should be distinguished 
from the clean case, where the DMFT approximation 
causes several problems such as the stability of the spin 
liquid state '36[ and strong dependence of the dimension 
of spin dynamics [Tt!]. 
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Appendix A: Derivation of Eq. ([3]) from Eq. ([1} in 
the Replica method 

The replica trick [s!] has been utilized for the disorder 
average, given by the following identity 

WZ = lim ^" ~ \ (Al) 

where O means the disorder average of an operator O. 
is the replicated partition function 



(A2) 

where the corresponding replica action is 

Q icra ijaa 

+ ^ + E,)dUr) + E Jn^^ ■ S ■ (r) 

iaa ija 

+ J2iv4:ir)dlir)+U.c.)] (A3) 



with the spin index a = 1, Af and the replica index 
a = 1, n. 

The disorder average for the replicated partition func- 
tion is straightforward, given by 



(A4) 



where P[ei\ is the Gaussian distribution function with 
J deiP[ei\ ~ 1. Performing integrals for random vari- 
ables, we obtain an effective action Eq. ([3]). 



Appendix B: Derivation of Eq. ^ from Eq. ([3]) in 
the cavity method 

We solve the replicated Anderson lattice model Eq. 
^ in the DMFT approximation. We apply the cavity 
method for the DMFT mapping [l6j 



(Bl) 



where S'^ is the part of the action at a particular site 0, 
ASn is the part of the action connecting the site with 
other sites, given by 



drJ2{''l»i9r - M)C(r) + dll{T)idr + i?,)dS.(r)} 



'2M 



dr / dr' ^ [cJ:(r)dgJr)-H4:(r)C(r)][ct(r')dL'(^')+4''(^')c[;.-(r')], 



(B2) 
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J 



icra 

2 rl^ 



2zM , 



(B3) 



respectively, and Sn '' is the rest part of the action. 
The partition function can be expanded as follows 



1- / dr(A£(T)) 





(0) 



/3 /-/S 



1 dr dr'(r.A£(r)A/:(r'))(") + ---(|?4) 



j drj dr' SS(r) • S^(r) Sg(r') • S^(r'), ^^^^^^ 

^ ^ iab 



with = /2?<^oa2'rf"^o.exp(-5i°^) and AS,, 

J.^drACir). 

The non-trivial term in the first order is given by 



(0) 



2zM 

7-2 /-/J /-/J 



^ f dr f dr'Y, SLpir)RfJ°s\r-r')S'o,5ir'), (B5) 



where <5)(r-r') = (T.^^Jr)^^,(r'))(°). The second 



order term is 



f^dr f dT'{TrAC{T)AL{T'))(') = [' dr /' dr' J] (r.4",(r)%C(r)4^(r')t,ocL(r'))^°^ 
Jo Jo Jo Jo I. 

= - f dr f dr' Y 4".(^)%GfJ°) (r - r%ocl{r'), (B6) 



■ijaah 



I 

where G'^y°^(r — r') = — (TVc^crl''")^]^ ('''')) One can Appendix C: Derivation of critical exponents, Ac, 
easily verify that all higher order expansions in Eq. (jB4[) 
vanish in the limit z — >■ oo, which is at the heart of the 
DMFT approximation [ll|. 



In the z — )■ (X) Bethe lattice we perform further 
simplification [l^ 



A/, and A(, 
1. At zero temperature 

Inserting Eqs. (gSj-dlSl) into Eqs. (I34l)-(l36l), we obtain 



V 



jab (0) , 



ah 



(B7) 



As a result, we reach an effective single-site action Eq. 
(g]) called the DMFT approximation. 



M 

2 /i3/-. . , ,3Af-l 



(CI) 



J^4C3A,-ic.^^^-\ (C2) 
Sb(c^) = V^A,A}Ca^+a,-i^^'+^'~\ (C3) 
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It is naturally expected 



2/(0+) = -Ed-X, 
Sfc(0+) = -A 



(C4) 
(C5) 



for power-law solutions at zero temperature. 

Combined with the Dyson equations ([37 |) - (p9l) . we 
reach the following equations 



Using the property of Ca-i = ACa, we obtain 

A/ = MAfc. 
As a result, the first solution is 

M 



A 

Afc - 



^ Af + 1 ' 
1 



Af + 1' 



(C13) 
(C14) 



A,Ca 



1 



M' 



Ae-l 



M' 



j24C3A,-l^'^^-\ 



(C6) 



(C7) 



-OJ 



-F^A,A/CA.+A,-lW^=+^^-^ (C8) 



ACa^-I 

The last equation gives 

A, + A/ + Ab = 2. (C9) 

Comparing this with the first equation, we get 

A, = 1. (CIO) 

The second equation gives two possible solutions. One 
is again Ac-I-A/-|-Af, = 2, and the other A/ = Ab = 1/2. 
We can find the first solution by equating the coefficients 
in Eqs. ((C7)) - (|C8)) and obtain 



exactly the same as those of the multi-channel Kondo 
effect. 

2. At finite temperatures 

Inserting Eqs. dSH]) -([521) into Eqs. we obtain 

Ee(iw/) = ^AfAtl3'-^f-^^^'fi,iiCdi) 

+ f-^p^-^^^,{^u^^l (C15) 

l^fiiuji) = ^A,Abl3^-^^-^'m,b{iooi) 

- J^A^p'-^^f^fffitLoi), (C16) 

Mm) = V^A,Af/3'-^^-^f^,f{ii?i), (C17) 

where 



1 



Caj-iCa^ 
1 



M' 



AcAfAb, 



-V^A,AfAb. 



These two equations result in 

C'Aj'-lC'At = A/CAi,-lCAj,. 



(Cll) 



(C12) 



dte'^^*gfit)gbi-t), (C18) 



(C19) 



dte^-^^[gfit)rgf{-t), (C20) 



^cf{m) = / dte"'^'gc{-t)gf{t). (C21) 
Jo 

Using the Dyson equations, we obtain the final self- 
consistency expressions 



+ j2^3^2-4A,[^^^^(^-^)_^^^^(^-^)]^ (C23) 

K^%\m) = inr^' - [A - S,(iPo)]/3'^^" + Y"" A,Af^''-^^-^^-^-\^,f(in) - *c/(«i^o)]- (C24) 



As the zero temperature case, we obtain two power-law A/ -f- Af, = 1, (C26) 

solutions, comparing the powers of /3 terms. One is ^.^j^ ^ ^^^^^ ^^^^^.^^ 

A, = 1, (C25) A, = 1, (C27) 
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A/=A6 = l/2. 



(C28) 



Note that Eqs. (|C22|) and (|C24|) are same for both sohi- 
tions. Only Eq. (jC23[) distinguishes these two solutions. 



Inserting Eq. ^ into Eqs. (Hg, (ICT8)) - (ICT9t . we 
obtain 



and 



^c{ii^i) = iTT sgn(a;i), 



r(i-A/) 



<^t{m) = (2^)^^(-i)' 



^fb{iuji) = ITT sgn(w;), 



r(i-^ + t)r(i-^-f)^ 
r(i-Afc) 

r(i-% + t)r(i-%-t) 



r(-A,) 



^cbitCoi) - (27r)^^+l^(-l)'^ ^ = ^ ^ ^, 

V2 2 ^ 27r/-^ V2 2 277-' 



^fffiiu^i) = (27r)3^/z(-l)' 



r(l-3A/) 



r(i-^ + f)r(i , 
r(-A/) 



3A/ _ (£j_N 

2;t' 



ri'i — =i -I- i^L'lFl'i — =J- — ) 
^2 2 2iT ' ^2 2 2itI 

I 



(C29) 
(C30) 

(C31) 



(C32) 
(C33) 

(C34) 

(C35) 



Inserting these expressions into Eq. (jC22p . we obtain and the equation 
the equation for Ac 



1 



From Eq. (IC24[) we obtain the condition 

1 



A$6(wo) 



A6N12 



[r(i - At) 

A(27r)Ai> (-l)"r(l- Afa) 



(C38) 



Inserting Eqs. ((C3T|) and ([C35|) into Eq. ((C38l) . we ob- 



(C37) tain 



rfi - Ai. + ^L^rd - ^ - 

^ V-^ 2 ' 2ir''-^V-^ 2 _ — 0) 



(-i)'r(i-Afc 

One can show that 



V^AcAfAb{2^) 



(-i)'+ir(-A;) 



V2 2 ^ 2;t'' V2 2 27r'' 



rn — Ai. -I- ]^L\v(^ — At _ ii."! 

^ 2 +2.MU 2 ^_(; = 0) 



(-imi-A,) 

(-i)'+ir(-A;) 



V 2 2 ' 2ir /-^ V 2 2 27r / 



r(i- Ab) 
r(-A/) 



nU(fc-H^) ^ 

T\k=ii^ ^ h 2^) 

" nU(fc-^ + ^) ^ 



]39) 



(C40) 
(C41) 



Here we have used A/ + Ab = 1. From Eqs. (IC39I) - (|C41|) we obtain 



1 = -V^AcAfAb{2TTf 



r(-A^)r(i-A,) 



r(i + ^)r(i-4t) 
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^V^A,AfAhr{-Af )r{l - Afc)[2sin(^A//2)]2. 

(C42) 



Equation (|C23p gives the condition 

1 



E/(iwo) - Ed - X = T 



l-Af 



r(|-^)r(i~^) 

A/(27r)^/ i(-l)T(l-A/) ■ 

(C43) 



and we will consider them separately. 

a. First solution: A/ + A;, — 1 with A/ > 1/2 
The scaling Eq. (|C23|) becomes 

A/[<i>f{iuJi) -$/(iwo)] = -—AcAb[^cbiii^i) - *c6(iwo)]. 



(C44) 



The scaling equation (|C23| distinguishes two solutions, Inserting Eqs. (|C30|) . (|C33| into Eq. (jC44|) . we obtain 



(-imi-A,) 

One can show that 



M 



A,AfAb{2nf 



-l)'r(-Afc) 



T(- — Ai _i_ i:iL\T(- — Ai. _ 

V2 2 2ir/-^ V2 2 2ir/ 



(/ = 0) 



r(i-^ + t)r(i^^-t) 
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From Eqs. (IC45]) - ((C47l ) we get 
1 . -^.4.A,A(2.)= r<-^')r(l-A/) 



A,AfAbT(-l^b)T{l - A/)[2sin(^A//2)]2 



(C48) 



Equations (|C42p and ()C48P result in the following 
equation 

Afr(-A/)r(i - Ab) - r(-Ab)r(i - a^), (C49) 

or equivalently 

AfA6 = A/. (C50) 
I 



I 

As a result, we obtain the solution 

M 

^ ~ A/ + 1 ' 
1 



Afc 



[I = 05(145) 

(C46) 
(C47) 



M + 1' 

with the condition A/ > 1, and M = 2 actually. 

b. Second solution: Af = At — 1/2 
In this case the scaling equation (|C23|) becomes 



(C51) 
(C52) 



V 



(C53) 



where both the hybridization and RKKY interactions Af, — 1/2, where the RKKY term of fff{iuji) will give 
give rise to the same order of magnitude for self-energy a similar result, 
corrections. We can use Eqs. (jC45p - (jC47p with A/ = 
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